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Abstract: Pair production in a constant electric field is closely analogous to bubble 
D . nucleation in a false vacuum. The classical trajectories of the pairs are Lorentz invariant, 



but it appears that this invariance should be broken by the nucleation process. Here, we 
use a model detector, consisting of other particles interacting with the pairs, to investigate 
how pair production is seen by different Lorentzian observers. We focus on the idealized 
situation where a constant external electric field is present for an infinitely long time, and we 
consider the in-vacuum state for a charged scalar field that describes the nucleating pairs. 
The in-vacuum is defined in terms of modes which are positive frequency in the remote past. 
Even though the construction uses a particular reference frame and a gauge where the vector 
potential is time dependent, we show explicitly that the resulting quantum state is Lorentz 
invariant. We then introduce a "detector" particle which interacts with the nucleated pairs, 
and show that all Lorentzian observers will see the particles and antiparticles nucleating 
at rest in the detector's rest frame. Similar conclusions are expected to apply to bubble 
nucleation in a sufficiently long lived vacuum. We also comment on certain unphysical 
aspects of the Lorentz invariant in-vacuum, associated with the fact that it contains an 
infinite density of particles. This can be easily remedied by considering Lorentz breaking 
initial conditions. 
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1 Introduction 

Vacuum decay through bubble nucleation was first analyzed in the ground-breaking paper 
by Voloshin, Kobzarev and Okun (VKO) [1]. Their analysis also uncovered an unexpected 
problem. In the semiclassical picture, the bubble is formed momentarily at rest; then it 
expands, rapidly approaching the speed of light. Lorentz invariance of the vacuum requires 
that the probability of nucleating a bubble should be the same in all inertial frames. This 
suggests that the rate of bubble formation per unit spacetime volume should include an 
integral over the rest frame of nucleation, that is, over the Lorentz group. This integral is 
of course divergent, so one obtains a meaningless infinite answer 1 . 

1 For a recent discussion of related issues, see [2-4]. 
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An elegant resolution of the paradox was given by Coleman [5], who developed an 
instanton method for calculating the bubble nucleation rate. The instanton in this case is 
0(4) invariant, and its Lorentzian continuation is invariant with respect to Lorentz boosts. 
From this Coleman concluded that "an expanding bubble looks the same to all Lorentz 
observers, and to integrate over the Lorentz group is to erroneously count the same final 
state many times" . 

There is little doubt that Coleman's calculation of the nucleation rate is correct. How- 
ever, some intriguing questions about the bubble nucleation process remain unanswered. 
The Lorentzian continuation of the instanton solution describes a bubble which contracts 
from infinite size, bounces at a minimum radius, and then re-expands. This solution is 
Lorentz invariant, but it appears that the contracting part of it is unphysical and needs 
to be cut off. However, the cutoff would necessarily break the Lorentz symmetry, and we 
would be back to the VKO problem. The answer that is usually given is that the semi- 
classical approximation breaks down near the nucleation point (that is, near the bounce at 
the minimal radius). So there is no well defined cutoff surface. The bubble emerges from 
a fuzzy quantum region, and only its late-time asymptotic behavior has classical meaning. 

This may be true, but we can still ask how different observers would see the nucleation 
process. An observer can fill space with detectors measuring the scalar field of the bubble. 
This would yield a complete spacetime history of the nucleation process. In particular, the 
rest frame of bubble nucleation may be determined in this way. 

One can anticipate several possible scenarios that can emerge from such an analysis. 
(A) It could be that the frame of nucleation simply coincides with the rest frame of the 
detectors. In other words, each observer will see bubbles forming at rest in her own rest 
frame. (B) It is conceivable that the contracting part of the bubble history is not completely 
cut off and can be at least partially observed. In fact, in the case of scenario (A), an observer 
moving relative to the detectors will see the detection begin on the contracting part of the 
bubble trajectory. (C) Yet another possibility is that the frame of nucleation is influenced 
by how the decaying false vacuum was set up. If the false vacuum has zero energy, its 
spacetime is flat, and the space will be filled by nucleating bubbles in a finite amount of 
time. This implies that the false vacuum could not have existed forever; it must have been 
created in some manner in the past. We could imagine, for example, that the false vacuum 
was set up on some spacelike surface, e.g., on a hyperplane t = const. The preferred 
frame of reference specified by this choice could in principle determine the frame of bubble 
nucleation. 

If the false vacuum has positive energy, it is described by de Sitter space, and the 
bubble nucleation process can continue forever. However, the geometry of the full de Sitter 
space is similar to that of the worldsheet of a bubble: it describes a universe contracting 
from infinite size, bouncing at a minimum radius and re-expanding. Bubbles would fill 
the universe during the contracting phase, so one needs to introduce a cutoff. This would 
break the de Sitter symmetry and define a preferred frame 2 . One might think that the 

2 The conclusion that the false vacuum state cannot be extended to the infinite past has been shown to 
hold in the most general case, without assuming homogeneity or isotropy [6, 7]. 
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memory of the preferred frame would gradually fade away and would have no effect on 
bubbles nucleating at late times. However, it was shown in [8] that a peculiar "persistence 
of memory" effect does not fade away and persists at arbitrarily late times. Thus the 
possibility (C) that the frame of nucleation is influenced by the initial conditions remains 
a viable option. 

In the present paper we shall attempt to address these issues using pair production 
in a constant electric field as a model for bubble nucleation. We shall consider charged 
spin-0 particles, described by a complex scalar field (j>(t,x), in a constant electric field E 
in a (1 + l)-dimensional Minkowski space. Such an electric field is Lorentz invariant, and 
pair production by the Schwinger process [9] is closely analogous to bubble nucleation. 
It is conceivable that despite the invariance of the background, the Lorentz symmetry 
could be spontaneously broken by the quantum state of the 0-field. Indeed, particles and 
antiparticles of the pairs are accelerated by the electric field in opposite directions, and 
since the number of the pairs is constantly growing, one can expect a growing expectation 
value of the electric current (J 1 ). Any non-zero value of (J M ) would break the Lorentz 
invariance. However, the situation is complicated by the fact that the current (J M ) is 
formally infinite and needs to be renormalized. 

Hence, the first part of our project is to investigate the invariance properties of the 
in-vacuum state, defined by requiring that the mode functions of the i^-field are positive- 
frequency at t — > — oo. This quantum state has been extensively discussed in the literature 
(see, e.g., [10-14] and references therein), but to our knowledge its invariance properties 
have not yet been studied. We show by a direct calculation that the in-vacuum state is 
Lorentz invariant. 

Nonetheless, as we shall see, divergent terms in the expectation values of physical 
observables, like the current or the energy- momentum tensor Ta V , do not seem to admit 
any Lorentz invariant regulators. The unusual properties of the in-vacuum state are due 
to the fact that at any finite time it contains an infinite number of out-particles. In a more 
realistic model, the electric field would be cancelled by the back-reaction of the created 
pairs, and only a finite density of particles would be produced. Alternatively, in the context 
where the electric field is external, we can still consider states where the number density of 
pairs is finite at any given moment of time. For instance, we could choose initial conditions 
such that the number density of pairs vanishes on some initial surface t = const. 3 We shall 
come back to this issue in Section 4. 

Despite its somewhat unphysical properties, we find that the in-vacuum is a useful 
laboratory for studying the bubble nucleation process, at least in the limit where the initial 
conditions are removed sufficiently far in the past. In this paper we shall consider tree-level 
interactions between the pairs and the detector. The kinematics of these interactions is 

3 The number density of pairs is not a sharply denned quantity. Nonetheless, it can be defined by using 
the instantaneous Hamiltonian diagonalization. Such states (with vanishing particle number at t = 0) have 
been discussed in Refs. [10] and [14]. The initial number of particles could also be rigorously defined by 
considering an electric field which vanishes at t — > — oo and then turned on at some time in the past. 
An abrupt turn-on at an instant of time was considered in Refs. [14-16], and an adiabatic turn-on with 
E(t) oc l/cosh 2 (t/T) was discussed in Ref. [10]. 
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such that, out of the infinite bath of created particles, only those whose momentum relative 
to the detector is in a certain range will have a chance to interact with it. In this sense, 
most particles in the bath are invisible, and their infinite density is irrelevant. 

To be specific, we model the detector by introducing two additional scalar fields: a 
charged scalar t[;(t,x) and a neutral scalar a;), with the interaction Lagrangian 

L iBt = -g{<f>Hx + ^<l>x), (1-1) 

where g is a coupling constant. This model has been studied in great detail by Massar and 
Parentani [12] and by Gabriel et al [13], who used it to investigate the Unruh effect for an 
accelerated detector. 

We start with a charged ^-particle in the initial state. It interacts with 0-antiparticles 
of the pairs via 4 ificj)* — > x an d thus has a finite lifetime rw,. We shall calculate the mo- 
mentum distribution of x-particles in the final state and use it to deduce the momentum 
distribution of the created </>-pairs. To achieve this goal, we will have to consider the inter- 
action (1.1) with a time-dependent coupling, g(t) = g exp(— t 2 /T 2 ), so that the detector is 
turned on for a finite period of time At ~ T. We shall also briefly discuss the case when 
the role of the detector is played by the neutral x-P ar ticle. All our results point in the 
direction of option (A) - that the frame of pair (and bubble) nucleation is determined by 
the frame of the detector. 

The paper is organized as follows. In Section 2, we introduce the in-vacuum state 
and review the Schwinger pair production using the method of Bogoliubov coefficients. In 
Section 3, we show that the in-vacuum is Lorentz invariant. In Section 4, we calculate the 
two-point function in this vacuum and discuss the expectation value of the current, pointing 
out the pathologies associated with the infinite density of particles in the in-vacuum. We 
argue that these can be remedied by considering Lorentz breaking initial conditions. In 
Section 5, we set up our detector model. The final distribution of x-particles and the 
observed momentum distribution of the pairs are calculated in Sections 6 and 7. Finally, 
our results are summarized and discussed in Section 8. Some technical details of the 
calculations are presented in the Appendices. 



2 Schwinger pair production 

We consider a constant electric field in (1 + l)-dimensions. Spin-zero charged particles that 
are being pair produced by the Schwinger effect are described by a complex scalar field 
4>(t,x); the corresponding action is 



S 



d 2 x 



-rT (9a + ieA^) <j>* {8 U - ieA u ) - m 2 0> - ^-F^ v F f 



1 



(2.1) 



Here, m and e are the mass and the electric charge of the particles, 



fits 



dn,A u d u A, 



4 Here and below we use asterisk to denote antiparticles. 



(2.2) 
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e^y is the unit anti-symmetric tensor with eoi = 1, and E = const is the value of the electric 
field. It is clear from this expression that a constant electric field in (1 + l)-dimensions is 
Lorentz invariant. More specifically, an observer with a 2-velocity sees the electric field 
given by 



E^l = Ffj, v u u = E , 



(2.3) 



where is a unit spacelike vector normal to u^. 

For this electric field, we may choose a Lorentz-covariant gauge in which the gauge 
field is given by 

A ii = X -Ee ilv x v . (2.4) 

However, since the calculation in this gauge seems technically more involved, we choose a 
non-covariant gauge where 



This gives the components of the gauge field as 

= (A t , A x ) = (0 , -Et) , 

where (t, x) = (x ,^ 1 ). 

2.1 The in-vacuum and the Bogoliubov coefficients 

The variation of the action with respect to (ft gives the field equation, 



(2.5) 



(2.6) 



-dl + (8 X - ieA a 



m 



0. 



We expand the field in terms of the creation and anihilation operators, 



(ft{t,x) 



dk 



(2vr) 1 /2 

where the mode functions cftkit) satisfy the equation 

~d 2 



a k <ft k (t)+bl k <ftl(t))e 



dt 2 



+ m z + (k + eEty 



0. 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



The canonical commutation relations lead to 

[a k , a\,} = 5(k - k') , [b k , b\,\ = 5(k - k') , 

and the normalization condition, 

i[<ftl(t)d t <ft k {t) - Mt)dt4>l(t)) = 1 • (2.11) 

Linearly independent solutions of Eq.(2.9) can be expressed in terms of the parabolic 
cylinder functions, 



{z)<xD v .[±(l-i)z], 



(2.12) 
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where 

r—fk\ 1 -HA m? , 

The general solution is given by a linear superposition of 4^. We choose 

MZ) = j2^ e ^ DA ~ {l ~ i)z] ' (2 ' 14) 

where the coefficient is chosen in order to satisfy the normalization condition (2.11). For 
future reference, we give a useful integral representation of the parabolic cylinder functions, 



e * 



,2 



D A (Z) = ~ I dwe^^-^w-^ 1 , for »(A) < . (2.15) 

r(— A) jo 

From the asymptotic expansion formula, 

z 2 3 
D p (z) ~ e~~z p , for |z|»l, \z\ 2> |p| , |argz|<-7r, (2.16) 

we find 

'V2\z\Y e^ 2 , for z < - \u\ . (2.17) 



(2e£)V4 

Hence idt4>k ~ —eEtipk at t — >■ — oo, indicating that these mode functions are positive 
frequency in the "in" region at t — > — oo. The corresponding vacuum state |0)i n , defined 
by afc|0)i n = 6fc|0); n = is the in- vacuum, which has no particles in the asymptotic "in" 
region. 

Because of the non-trivial background, the positive frequency mode function <f>k at 
t —> — oo does not remain positive frequency at finite t, and in particular it is given by 
a linear combination of the positive and negative frequency functions at t — >■ +oo. The 
positive frequency functions at t — > ±oo are related by a Bogoliubov transformation, 

^ = a fc 0r + /« Ut % (2-18) 

where 

^ Ut(z) = pjv? e " lf " + l)z] ' (2 ' 19) 

is the positive frequency mode function at t — > oo and \ak\ 2 — |/?fc| 2 = 1- We can check that 
the asymptotic expansion of 0£ ut indeed gives 

1 

' k ~ (2eE)V* 



f y/2z) e"5 2 , for * » . (2.20) 



Using a linear relation, 

ZV[-(1 - *)*] = e^*A,*[(l - i)z] + f7§) e " ! ^[( 1 +^ > ( 2 - 21 ^ 
we can read off the Bogoliubov coefficients, 

y/2^ 



Oik 



^~ v ) , p k = e™ v * . (2.22) 



r(-i/*) 

This gives 

|/3 fc | 2 = e-^ = e-^#. (2.23) 
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2.2 Momentum distribution and Lorentz invariance 

The Bogoliubov coefficients fik are simply related to the momentum distribution of particles 
in the "out" state at t — > oo, 

dn 1 , „ ,o , 

* = . (2 ' 24) 

This can in turn be related to the rate of the Schwinger process, ie, the rate of pair creation 
of charged particles. The mixing between positive and negative frequency modes in the 
in- vacuum mode functions (2.14) occurs in the interval 

M<M, (2.25) 

centered at z = (or k = —eEt). Hence, the number density of particles created between 
time to and t > to is given by 

i r—eEto 

n = — dk\h\ 2 . (2.26) 

2tt J_ eEt 

This leads to 



dn ( eE 
~dt ~ V 2vr 



(2.27) 



This is the Schwinger formula for the rate of pair creation. 

An important feature of the momentum distribution (2.24) is that it is independent of 
k. The total density of created particles, obtained by integration over k, is infinite. This 
is not surprising, since the pair creation process was going on at a constant rate for an 
infinite time. A more realistic calculation would include back-reaction of the pairs on the 
electric field, so the field would gradually decrease and only a finite density of pairs would 
be produced. But according to Eq. (2.27), in a weak electric field with \eE\ <C m 2 pair 
creation is a very slow process, and the field can remain nearly constant for a very long 
time. Thus, one can expect that our idealized treatment should apply in some limiting 
sense. We shall see, however, that taking the limit and interpreting the result is not always 
straightforward. 

It should be noted that k in Eq. (2.24) is the canonical momentum, which is related 
to the physical momentum by 

fcphys = k + eEt - (2-28) 
The only subtlety is in determining the range of this distribution. Pair creation in a given 
mode occurs in a spacetime region where the negative frequency contribution to the mode 
function becomes significant. As we discussed, this occurs at k ~ —eEt. Hence, Eq. (2.24) 
is valid in the range 

- eEt < k < oo. (2.29) 

The uncertainty in the lower limit of this range can be estimated from (2.25), which (as- 
suming A = m 2 1 (eE) 3> 1) leads to the uncertainty 

Ak min ~ A 1/2 m , for A»l. (2.30) 

Note that at future infinity, the distribution of out particles is given by (2.24) in the full 
range — oo < k < oo. The final distribution is then Lorentz invariant 5 . On the other hand, 

5 This conclusions follows from the Lorentz invariance of the phase space element dkdx [17]. 
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at any finite time, the invariance is broken by the lower cut-off in Eq. (2.29). 
Eq. (2.28) leads to d/c p h ys = dk, and since (3^ are independent of k, we have 

^ = f'*' 2 ' < 2 - 31 > 

From (2.28) and (2.29), the range of physical momentum for the created particles is 

< fcp h ys < °o. (2.32) 
Hence, we expect to see particles which predominantly have positive momentum 

^=4lA| 2 «(W), (2-33) 

and similarly 

^\Pk\ 2 9(-k phys ) (2.34) 



dk phys 2ir 

for antiparticles. This restriction on fc p h y s is intuitively plausible. Pairs are produced at 
rest; then particles and antiparticles are accelerated by the field in opposite directions. In 
the limit of infinite time, one can expect to find particles having arbitrarily large momen- 
tum, but one would not find particles with k p h ys < (or antiparticles with fc p h ys > 0) 6 . 

A potential problem with this picture is that the cutoff at /c p h ys ~ clearly breaks 
Lorentz invariance. This is not necessarily impossible. On one hand, even though the 
background is Lorentz invariant, the invariance could be broken by the in-vacuum state 
which we are using for the <^>-field. This would correspond to option (C) in the Introduction, 
where the frame of nucleation might be influenced by how the decaying vacuum is set up. 
Finally, we must consider the possibility that the in-vacuum is Lorentz invariant, but all 
observers see the distributions (2.33) and (2.34) for the momenta of particles and anti- 
particles relative to the detector's rest frame. This would correspond to option (A) in the 
Introduction. 

Our next task is therefore to investigate the Lorentz invariance of the adiabatic in- 
vacuum. 

3 Lorentz invariance of the in-vacuum 

To examine the behavior of the in-vacuum under Lorentz transformations, we compare two 
in-vacua defined with respect to different rest frames which are related by a Lorentz boost. 
To do so, we consider a gauge transformation that relates the components of the gauge 
field in the two frames, such that has the same form given by Eq. (2.6) in each frame. 
A positive frequency mode function of (p at t — > — oo in the gauge (2.6) reads 

Mz)e lkx ; z = vWi + ^V (3.1) 



6 Note, however, that the uncertainty in the lower limit for fc p hy S is of order Afc p h ys ,min ~ A 1/,2 m. For 
A > 1, this seems to leave some room for highly relativistic particles with both signs of the physical 
momentum. 
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with 4>k(z) from Eq. (2.14). We consider another Lorentz frame given by 

i=j(t-vx), x = 7O -vt); 7 = (1 - v 2 )" 1/2 . (3.2) 
In the barred frame, if we adopt the gauge Aq = 0, the positive frequency modes at i — > — oo 



are 

Akx 



where the prefix '0' means the gauge 

A- = 0. (3.4) 
Note that the functions o4>k{z) are the same as (j>k(z) with the replacement z — > z , 

ofoC*) = ■ (3.5) 

To compare the positive frequency function in the barred frame with the one in the 
original unbarred frame, we have to express the barred frame positive frequency fuction in 
terms of the original coordinates (t,x). To do so, we perform a gauge trasformation from 
Aq = to Ao = gauge. 

The electromagnetic vector potential trasforms from the barred frame to the original 
frame as 



dx u 

o\ = —A* , (3-6) 



dxv- 



where qA^ is the components of the potential in the gauge (3.4) in the coordinates (t,x). 
This gives 

o^o = ^rM = ivEi= 1 2 vE(t - vx) , 
at 

dx — 

Q A l = —A- l = - 1 Ei=- 1 2 E{t-vx). (3.7) 
We perform the gauge transformation from this gauge to the gauge (2.6), 

Afj, = A^ - d^n , (3.8) 
so that An is given by (2.6). Namely we set 

A = A - d tt = 0, Ai = A l -d 1 n = -Et. (3.9) 
This can be easily solved. We find 

J) = (t 2 - 2vtx + x 2 ) . (3.10) 

Under this gauge transformation, since the covariant derivative of </>, 

(0 M -ieA M )0, (3.11) 
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is gauge-invariant, the mode function trasforms as 

ofa^fa = ofae- ie(l . (3.12) 

Thus the positive frequency mode function in the barred frame is expressed in terms of the 
original coordinates as 

^e^; 4>- k (t,x) = Mz)e- ieQ(t ' X \ (3-13) 

where 

J. 

z = veE^{t — vx) H — -j= , x = j(x — vt) . (3-14) 

A simple way to examine the Lorentz invariance of the in-vacuum is to evaluate the 
Klein-Gordon inner product 

K= ! dx{<p k d t h - 4>- k dtfa ) e ~ ikx+fks (3.15) 

J f =const 

This product characterizes the mixing of positive with negative frequency modes. If K is 
non- vanishing, it means that the vacua defined by <p k and by 4>t are intrinsically different, 
hence the in-vacuum is not Lorentz invariant. The inner product (3.15) is computed in 
Appendix A using the exact mode functions, and the result is 

K = 0. (3.16) 

This shows that the in-vacuum (hence also the out- vacuum) is Lorentz invariant. 

We shall see in the next Section that this state has certain unphysical features. These 
would actually make it pathological when we consider interactions beyond the tree level. 



4 Two-point function, Current, and Lorentz breaking quantum states 

A quantity of particular interest for our problem is the expectation value of the electric 
current (■/«)■ The physical picture of particles and antiparticles produced at a constant 
rate and driven in opposite directions by the electric field, suggests that there should be a 
steadily growing electric current, 

^-C (4,, 
where C = const, while the charge density remains equal to zero, 

(J ) = 0. (4.2) 

On the other hand, the expectation value of any vector in a Lorentz invariant quantum state 
should be zero. Since we found that the in-vacuum is Lorentz invariant, this indicates that 
we should have (</«) = 0. Once again, there is some tension between Lorentz invariance 
and the physical picture of pair production. Let us investigate this question in some detail. 
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Formally, the current can be expressed as 

J M = | (tfDpt - 4> (D^) + h.c. (4.3) 

A naive calculation of the expectation value of in the in- vacuum gives an infinite answer. 
One way to regulate this sort of infinity, which respects Lorentz and gauge invariance, is to 
use the point splitting method. The idea is to take the two field operators in the products 
in Eq. (4.3) at different spacetime points, x and y, and then take the limit y — > x. The 
current expectation value ( J At ) can then be expressed in terms of the 2-point function, 

(?W = G + + G~ , (4.4) 

where 

G+(xV) = in (0|<£V) e -feJ?^" 0(^)|O) in , (4.5) 

G-(^,^) = in (O|0(^) e -* e ^ A " dxV </>V)|0) in , (4.6) 

and the Wilson line has been inserted to make the 2-point function gauge invariant. The 
expectation value of the current is given by 

lj ) = — h m — \G w (y y - x u ) = ie hm -A -, (4.7) 

where Ax u = y u — x u . 

The 2-point functions G iz {Ax^ 1 ) are calculated in Appendix B. The result is 

lal 2 W iX/2 (ieE As 2 /2) 

G+(Axn = ^T{- V *)^p===L± ) for Ai-Az>0 (4.8) 
4vr y/zeEAs 2 /2 

and 

lal 2 W_ iA/2 n (-^As 2 /2) 

G+(A^) = ^-r(-i/) A ° l =-^-, for At-Ax<0. (4.9) 

V ' 4vr v ' ^-ieEAs 2 /2 K ' 

Here, W a> p (z) are Whittaker functions, As 2 = —(At) 2 + (Ax) 2 , and a is the Bogoliubov 
coefficient given by Eq. (2.22). We have omitted the index of momentum k, since a is 
independent of k. The function G~ is simply related to G + , 

G~(At, Ax) = G + (-At, Ax). (4.10) 

For At + Ax < this is given by (B.9), and for At + Ax > it is given by (B.10). 

The first thing we note is that the 2-point functions G^ and are Lorentz invariant. 
This is a further manifestation of the Lorentz invariance of the in- vacuum. The function 
G^ is also invariant with respect to time reversal, 

G (1) (Af,Ax) = G (1) (-At,Ax), (4.11) 

but not with respect to spatial reflection, 

G (1) (Ai,-Ax) = G (1) (At,Ax)*. (4.12) 
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The latter is not surprising, since the electric field defines a preferred direction on the 
x-axis. G*- 1 ) is invariant under a simultaneous transformation Ax — > —Ax, E — > —E. Note 
also that the functions are not invariant under time reversal; hence the in-vacuum state 
is not time reversal invariant, as one might expect. 

For the calculation of the current (4.7), we are interested in the 2-point functions in 
the limit of Ax^ — > 0. This is given by 



G+(Ax^ 



47T 



, ieEAs 2 \ ,/l iX . 



(4.13) 



for At — Ax > and by a complex conjugate expression for At — Ax < 0. Here, 7 is Euler's 
constant. The leading term in the limit of large mass, m 2 3> \eE\, is 

^(A^-M^n (414) 

The singularity structure of the 2-point functions G^, and therefore of G^ is rather 
unusual. Most importantly, the leading divergent term at As 2 — > is different from that 
for a free field: the coefficient multiplying the logarithm in (4.13) has an extra factor of 
|a| 2 , compared to the case of E = 0. In terms of the standard terminology, the 2-point 
function is not of the Hadamard form. Another non-Hadamard feature is the ^-function 
discontinuity arising from the different forms of G ± at At = ±Ax. 

The non-standard singularity structure of the 2-point functions makes the regulariza- 
tion of the expectation values of the current and of the energy-momentum tensor T^ w 
problematic. The point splitting method with Hadamard subtraction cannot be used, and 
even if one does not insist on the Hadamard form of the 2-point functions, the coincidence 
limits of G^ and of its derivatives are ill defined, because of the discontinuity on the 
light cone. The Pauli-Villars regularization also appears to fail, since the coefficient mul- 
tiplying the leading divergence in (4.13) is mass-dependent. This seems to indicate that 
Lorentz invariant regularization of the expectation values of observables is impossible in 
our in-vacuum. 

As mentioned in the introduction, this is to be expected, since the Lorentz invariant 
in-state contains an infinite number of created pairs. In other words, the divergent current 
is caused by "actual" particles, corresponding to the non-vanishing occupation numbers of 
all momentum states with k + eEt > (rather than by zero point fluctuations, whose effect 
can be subtracted by using Lorentz invariant counterterms) . 

In Appendix C we calculate the expectation value of using a direct momentum 
space regularization, 

\k\<.k max , (4.15) 

and then taking the limit A: max —> 00. Physically, the current is due to the created pairs, 
and therefore the regulator (4.15) amounts to ignoring the contribution from pairs with 
\k\ > k max . Noting that the mixing between positive and negative frequency modes occurs 
near the time when fc p h ys = k + eEt = 0, we can think of (4.15) as a proxy for the situation 
where the electric field is only turned on for a finite interval of time, \t\ < fc max /(ei?). The 
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regulator (4.15) is not manifestly Lorentz invariant and, not surprisingly, the result which 
is obtained by using it is not Lorentz invariant either; it agrees with Eqs. (4.1), (4. 2) with 

„ e 2 E T m 2 

C= e~— . (4.16) 

7T 

This is what one would expect on physical grounds: the rate of current growth is propor- 
tional to the rate of Schwinger pair production. 

In view of the unphysical properties of the Lorentz invariant in-vacuum, we are led to 
the conclusion that (in the presence of a constant electric field), no physical state can be 
Lorentz invariant. 

To illustrate this point more explicitly, we may consider the closely related example of 
massless spinor QED in 1 + 1 dimensions [18]. Aside from the conserved electric current 
= etp^ifj, this model also has an axial current, which in 1 + 1 dimensions is related to 
the electric current by 

jn = e ^ 7 M 7 5^ = ^u Jv _ ( 417 ) 

Here, 7 5 = 7°7 1 . As is well known the axial current would be classically conserved, but 
any gauge invariant (and Lorentz invariant) regulator will lead to a nonzero divergence 
given by the anomaly equation [19, 20]: 

e 2 E 

V" = . (4.18) 

7T 

This equation is of course valid even if the electric field is not constant, but here we are 
primarily interested in the case where it is, and hence the background is Lorentz invariant. 
Combining Eqs. (4.17) and (4.18), the renormalized electric current satisfies 

e 2 E 

e^J„ = — . (4.19) 

7T 

This is consistent with our Eqs. (4.1), (4. 2), (4. 16) in the massless limit m = 0. Note that 
even though Eq. (4.19) for the current has a covariant form, it does not admit Lorentz- 
invariant solutions. In this sense, Lorentz invariance is spontaneously broken: quantum 
states satisfying Eq. (4.19) cannot be Lorentz invariant. 



5 A model detector 

We now turn to the question of how the pair production process is seen by different ob- 
servers. To this end, we introduce a model detector which involves two additional scalar 
fields: a charged scalar tp(t,x) and a neutral scalar x(t, x ), with the interaction 
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d 2 x 



4>Hx + ^<f>x 



dtH 



mt ; 



(5.1) 



where g is a coupling constant. We denote the masses of (f>,ip and % particles by m^, and 
m x , respectively. The charge of ^-particles is the same as that of 0-particles, as required 
by gauge invariance of the action. 
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Suppose first that we have a charged ^-particle in the initial state. It will scatter on a 
0-antiparticle of the pairs, producing a x in the final state. Such an interaction represents 
a "detection" of </>*, with <fi* playing the role of a domain wall in our model. Alternatively, 
we can have a neutral x particle in the initial state, in which case the detection process 
yields a i/'-particle in the final state. 

It will be useful to consider the kinematics of the scattering process. Starting with 

W -> X, (5-2) 

let q, k and p be the momenta of ip, <j)* and Xt respectively. Momentum conservation 
requires that 

q + k = p. (5.3) 

We remind the reader that k is the canonical momentum of 4>, which is related to the 
physical momentum by feph ys = k + eEt. Similarly, for ?/>-antiparticle, q p h_ y s = Q — eEt. For 
the neutral x-particle, p p h ys = P- Hence, Eq. (5.3) implies 

9phys ~i~ ^phys — Pphys j (5-4) 

which simply states the momentum conservation at the moment of collision. 
In the absence of an electric field, we also have energy conservation, 

\J 9phys + m l + \j fc phys + m l = \JP 2 + m l- ( 5 - 5 ) 

Squaring this equation and combining with (5.3), we obtain, after some algebra, 

2m% k phys = M 2 q phys ± y/ql hya + ml ^M±-±m\m\, (5.6) 

where 

M 2 = ra\-m%- m\. (5.7) 

For a given momentum of the probe g p h y s , there are two values of /c p h y s for which scattering 
is possible: one corresponding to the (^-particle coming from the left and the other from 
the right. Note also that the scattering is kinematically allowed only if 

m x > m^p + m<£. (5.8) 

Similar relations can be written for the neutral probe scattering X'P ~^ V 7 - The correspond- 
ing relation for the masses is > m x + m^. 

In the presence of an electric field, the energy conservation (5.5) does not strictly apply. 
One way to think about this is that the interaction process is not momentary and involves 
some time uncertainty 5t. Charged particles are accelerated by the electric field during 
this time interval, so their energy is changed. Alternatively, energy non-conservation can 
be understood in terms of the Unruh effect [12, 13]. The charged detector particle ip moves 
with an acceleration a = eE/m^p and is exposed to an effective temperature T = a/(2ir). 
Processes with energy fluctuation AE can therefore occur with a Boltzmann suppressed 
probability oc exp(AE/T). For example, a V _ P ar ticle can decay, ip — > (fix, even when 
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(a) (b) (c) 

Figure 1. Interaction processes in the <fi4>x model with a ^-particle in the initial state, (a) Decay 
of tp into <p and In the presence of an electric field, this can occur even if m^p < + m x . (b) 
True detection process: if) collides with a <fi* coming from a </></>* pair, producing a x- ( c ) A triplet 
of ijj* , <j) and x-particles is spontaneously produced from the vacuum. The initial -0-particle does 
not participate in the interaction. 

rrhip < nicf, + m x , so the decay would be kinematically forbidden in the absence of an 
electric field. 

Apart from the 0-pair production, ^-particle pairs will also be produced by the electric 
field. In addition, triplets of 4>il>*X an d of 4>*ipx wm De spontaneously produced from the 
vacuum. Various interaction processes in our model are illustrated in Fig. 1. In order to 
avoid false detections, we have to choose the parameters of the model so that processes 
shown in Figs. 1(a) and 1(c) are suppressed compared to the true detection process shown 
in Fig. 1(b). This can be achieved by requiring that 

m^p » m,0 ; (m x — m^)m^ » . (5.9) 

The first and the second of these conditions ensure the suppression of processes 1(a) and 
1(c), respectively. 

We shall calculate the amplitudes for the detection processes like (5.2) to lowest order 
of perturbation theory in the coupling g, but exactly with respect to the background electric 
field E. We first expand all the fields in terms of the in-vacuum creation and annihilation 
operators as we did in (2.8), 
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Here, </>&, ip q and Xp are positive frequencty mode functions of the respective fields at 
t = — oo. For (f> and tp fields they are given by Eq. (2.14) with masses and my,, 
respectively. For the x held, 

Xp = ^==jP*-^\ (5.11) 



where uj p = yp + m|. Since % is neutral, Xp 1S a ^ so positive frequency at t = +00. The 
canonical commutation relations are 

[a k ,al] = 5(k-k'), [b k ,b{,] = 5(k-k'), 
K, 4,] =%-</), IU fl<) =S(q-q'), 
[cp,cl,} = 5{p-p'). (5.12) 

In the interaction picture, the initial state | i) at t — > —00 evolves to a final state | f ) 
at t — > +00, which is given by 

/oo 
dt'H int (t')} |i> 
-00 

/■oo 

~ I i> — * / cft'flintCOli), (5.13) 



where in the second line we kept only terms up to linear order in g. 
With a ■(/'-particle detector, the initial state is 

|g)=4|0), (5.14) 

where | ) stands for the in- vacuum state, and we can examine the detection process by 
studying the distribution of x-particles in the final state, 



dN x _ (f I cpc p \ f) 



dp (q\q) 

To the lowest order in g, the out-expectation value in (5.15) can be expressed as 



(5.15) 



dt' I dt" / dki d Pl / dk[ dpi 

-OO J— OO J J 



x (q-lflint^lpiX) {p[,k[\ d p Cp \pi,h) (pi, ki\H int (t')\ q ), (5.16) 
where we have inserted two sets of intermediate states defined in the in-vacuum, 

ip',fc , ) = 4'4'i°)- ( 5 - 17 ) 

Strictly speaking, we also had to include intermediate states of the form 

4 4 c l' a i> 1 ) ( 5 - 18 ) 

and 
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which would also give non- vanishing matrix elements of Hi Qt . Inclusion of these states 
would account for the 4>ip*x an d 4>*i>X triplet creation processes of Fig. 1(c). We disregard 
these processes, assuming that the conditions (5.9) are satisfied. 
Using Eqs. (5.12), we have 



(p'l, k[\ 4> c p IP1> = 6 & ~ Pi) S (P 1 ~ P) 6 ( kl ~ 
and Eq. (5.16) simplifies to 

(f| c\c v |f) = J dk \M(p,q,k)\ 2 , 

where 



M(p,q,k) 



dt' {p,k\H int (t')\q) 



Substituting the mode expansions (5.10) into Eq. (5.1) for H{ n t, we obtain 



M(p,q,k) 
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2vr 



5(q-p-k) / dt' (plt/jqxl ■ 



(5.20) 



(5.21) 



(5.22) 



(5.23) 



Here, the delta-function enforces the momentum conservation, p = q + k. Substituting this 
in Eq. (5.21), we find 
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5(0) 



dt 4>* q _ p tpq X*p 



( f | c p c p I f ) 

Finally, substituting into (5.15) and taking into account the normalization 

(q\q) = 5(0) , 

we obtain the momentum distribution of ^-particles in the out-state, 



dNx 
dp 



1 

2^ 



dt 4>*_ p tp q Xp 



1 

= 2^\A x (p;q)\ ■ 



(5.24) 



(5.25) 



(5.26) 



In the alternative setup, with a neutral ^-particle in the initial state, the final distri- 
bution of -0-particles is given by 



dNj 
dq 



(f | d° utt d° ut 
(p\p) 



f: 



(5.27) 



where \p) = c P \0) and d° utT , d° ut are the creation and annihilation operators for ip- 
particles in the out-vacuum, which are related to the in-vacuum operators by Bogoliubov 
transformations. Assuming that 3> m^, ip-paii production can be neglected, so we can 
disregard the distinction between the in- and out-vacua for tp. Then, following the same 
steps as above, we obtain 

2 



tt 



dN„ 



dq 



1 

2^ 



9 / dt 4> q _ p ip q x P 



1 
2^ 



\Al){T,P)\' 



(5.28) 
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6 Detector lifetime 



Integrals of the mode function products A x (p;q) and A 1 p(q;p) appearing in Eqs. (5.26) 
and (5.28) have the meaning of amplitudes for the detection processes ip(fr* — > x and 
x4> ip, respectively. These integrals can be evaluated using the method developed in 
Refs. [12, 13]. Relegating details of the calculation to Appendices D and E, here we only- 
present the results. 

For a ^-particle detector, we are interested in the ^-particle distribution in the final 
state, given by Eq. (5.26). The amplitude A x (p;q) appearing in that equation is given by 



\A x (p;q)\ = — 



7T 



1/2 



fj\ 



V2eE 



A, 



(6.1) 



Here, uj p = ^Jp 2 + m 2 , is the x-particle energy, 
coefficient for (/>-particles, 



is the negative-frequency Bogoliubov 



exp 



eE 



(6.2) 



and \a^\ ~ 1 is the positive- frequency Bogoliubov coefficient for '(/'-particles. Furthermore, 



A^, =t A,, 



a± = 



(6.3) 



2 

mf 



(6.4) 



with i = (ft,ip,x, and W a>p (z) is the Whittaker function. Note that the amplitude (D.8) is 
independent of the ^-particle momentum q. 

This expression has a simple asymptotic form in the limit of large m x , 



A x S> A^A^. 



In this case, 



P.2 



2CT+, IU — 



I: 



(6.5) 



(6.6) 



where we have used an asymptotic representation W\^{z) ~ e Z//2 z A , which applies for 



\z\ > 1/T 



A 2 |. The x-particle distribution (5.26) is then given by 

.2 



dN- 



<r 



V 



4eEm 2 

X. 



exp 



eE 



dp 



(6.7) 



Note that the coupling g in our model has dimension of mass squared, so the right-hand 
side of (6.7) is dimensionless, as it should be. 

If in addition to (6.5) we assume that 3> m^, then the false detection processes of 
Figs. 1(a) and 1(c) are suppressed, and the x-particles in the final state count predominantly 
true detections. As we shall see in a moment, the detection rate is proportional to dN x /dp. 
Eq. (6.7) shows that this rate is proportional to the pair production rate of ^-particles, as 
expected. 
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In the general case, the amplitude A x (p; q) can be represented as 

_i 

I -A x (p; q) | = uj p 2 C x ( m x , m^, ) . (6.8) 
The corresponding ^-particle distribution is 



dN = B£^. (6 .9) 

The factor dp/ui p on the right hand side can be recognized as the volume element of the 
Lorentz group. This indicates that the distribution (6.9) is invariant under Lorentz boosts. 
The final distribution (6.9) can be represented as 

dN x f , dN x . . 

- ' ^ , (6.10) 



dp J dpdr ' 

where dN x /dpdr is the expectation number of x-particles emitted in the momentum inter- 
val dp per unit proper time r along the detector trajectory. The latter quantity transforms 
under Lorentz boosts as uj~ ; hence we can write 



dp dr uj p 

where the function / is Lorentz invariant and can therefore depend only on the invariant 

/ = e^p^ hys (r). (6.12) 

(The other invariant, I' = P/i9phys( T ) ^ s re l & ted to I by I' 2 = I 2 + m^rn^.) The physical 
momentum q p hys{T~) specifies the rest frame of the detector at the time of measurement 7 . 
Note that in this frame the invariant I is simply proportional to the x-particle momentum, 
I = -m^p. 

For a hyperbolic detector trajectory, we have 

dr = -^, (6.13) 

0>ip ^q,phys 



where = eE/m^ is the proper acceleration and ^ 9lP hys = y 9ph ys + m ^ * s * ne detector 
energy. Substituting this in (6.10) and comparing with (6.9), we can identify 

Cxi 2 _ 1 f dq phys /(/) _ (6 M) 



2VT J Wq,phys 

On the other hand, it follows from Eq. (6.11) that 



dNx 
dr 



f-f(I)- (6.15) 
J Up 



7 Here we have assumed that the distribution (6.11) depends only on the local value of the detector 
momentum q p h yB (T)- This is justified, as long as q p h yB remains essentially unchanged during the time of 
interaction. 
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and it is easily understood that the integrals over the Lorentz group in the last two equations 
are equal to one another. Hence, the scattering rate in the detector's frame can be expressed 
as 

^ = (6 . 16) 

(It 2tt m^p 
The characteristic collision time is 

^ = \C^E- (6 - 17) 

This is the average proper time that it takes for a detector to experience a collision. It can 
be called the detector lifetime. As one could expect, this lifetime is independent of the rest 
frame of the detector: it is the same for all points on the hyperbolic detector trajectory. 

For a x-particle probe, the detection process is x4> ~^ ip- The amplitude appearing 
in Eq. (5.28) is evaluated in Appendix E. The result is rather cumbersome, and we do not 
reproduce it here. As before, the amplitude is independent of q and can be represented as 

_ 1. 

\-A^(q;p) | = uj p 2 C^(m x , m^, m^) . (6.18) 
The corresponding detector lifetime is 



2irm x 
\a,,\ 2 eE 



(6.19) 

This is the same as 7V, in Eq. (6.17) with the replacement mwj — > m x and C x —> C^. Once 
again, we see that the detection rate is independent of the detector's rest frame. 

7 Momentum distribution of the pairs 

Our results in Section 6 indicate that the detection rate of </></>*-pairs is independent of the 
state of motion of the detector. This is in agreement with the Lorentz invariance of the 
in-vacuum state that we established in Sec. 3. 

However, the analysis of Sec. 6 does not allow us to determine the momentum distri- 
bution of the pairs measured by the detector. Indeed, it is clear from that analysis that it 
can only determine the distribution integrated over momenta or (which is essentially the 
same) over the proper time along the detector trajectory. In particular, we do not have 
an answer to the question we asked in the introduction: are the (j) and <jf particles of the 
pairs observed to nucleate at rest and then move away from one another, or they can also 
be detected on the 'wrong' parts of their hyperbolic trajectories, where they move towards 
one another at a high speed? In the former case, the distribution should be cut off at 
fcphys ~ 0, as in Eq. (2.33). 

In order to address this question, we consider a detector with a time-dependent coupling 

g{t)=ge- i2 ' T \ (7.1) 

so the detector is effectively turned on only for a finite time interval At ~ T around t = 0. 
Focusing on a charged ^-detector, we shall assume that 

<C eET <C , (7.2) 
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so that the detector momentum £/ p hy S remains essentially unchanged during this interval 
(while ^-particles can be significantly accelerated). For example, if we choose q = 0, then 
9phys ~ during the interaction. The resulting x-particle will then have (approximately) 
the same momentum as the </>*-particle that collided with ip, 



V ~ ^phys • 

We assume as before that 

In addition, we may assume that turning on and off of the detector is adiabatic, 



(7.3) 
(7.4) 

(7.5) 



so that the time variation of the coupling does not lead to a sizable violation of energy 
conservation at the moment of collision. In fact, the first strong inequality in (7.2) implies 



(7.6) 



so (7.5) is automatically satisfied in the regime where pair nucleation is rare, = m^/ (eE) 3> 
1, in which we are now focusing. It then also follows from Eqs. (7.2) and (7.6) that 

eET 2 > A > 1 . (7.7) 

To find the distribution of x~P ar ticles, we need to calculate the amplitude 

/oo 
dtg(t)<t>%(t)Mt)x* P (t). (7.8) 
-oo 

Since the ■(/'-particle is practically at rest during the time interval of interest, we can use 

1 



-IJTUi 



v / 2m^ 



(7.9) 



The amplitude (7.8) is calculated in Appendix F. Here we are only interested in the p- 
dependence of the amplitude. In the parameter regime specified above, this is given by the 
factor 



dp 



oc |.4 x (p; q = 0)| 2 tx exp 



(p + u p - m^,) 
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{eET? 



This distribution is peaked at p = p, which is specified by 



oi- 



and has width 



p + \/p 2 + m| = 



Ap ~ eET . 



(7.10) 

(7.11) 
(7.12) 
(7.13) 
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Note that we will only see a sharp peak in the distribution provided that \p\ 3> Ap ~ eET. 
Using (7.2), this requires 

f->m\. (7.14) 

Eq. (7.12) is to be compared with Eq. (5.6) for the kinematically allowed values of p. For 
(7 = 0, it gives 

pW± 2^ (m ' " m ' ) - (7 ' 15) 

In deriving (7.15), we have assumed — rr?^ 3> m^m^, which in turn implies the ultra- 
relativistic motion of the 4> particle, as in (7.14). 

From (7.12), it is clear that only one of the two possible values in (7.15) is actually 
detected. This is consistent with the picture of pairs being produced at rest (in the ob- 
server's frame) and then being accelerated by the electric field. 0-antiparticles move in 
the direction opposite to E and hit the V'-probe with fe p hy S < 0. This is imprinted in the 
momentum of the final ^-particle, p ~ A^ys < 0. 

From (7.14), we also conclude that the asymmetry in the momentum distribution is 
only visible for relativistic eft particles, which means that the precision to which we can 
determine the frame of nucleation is limited by the acceleration time, 

At nuc ~^. (7.16) 
eE 

This scale is comparable to the size of the instanton. 



8 Conclusions 

Our results indicate that all inertial observers in a (1 + l)-dimensional Minkowski universe 
with a constant electric field detect the same distribution (2.33), (2.34) of charged particle 
pairs created by the Schwinger process. In each observer's frame, particles and antiparticles 
nucleate at rest and are then accelerated in opposite directions by the electric field. This 
implies that particles that have already been formed and are ready to interact with a 
particular detector, may not yet exist from the point of view of another detector which 
moves relative to the first 8 . This picture fits well with the Lorentz invariance of the in- 
vacuum state which we established in Section 3. 

The situation here is similar to that for the Bunch-Davis vacuum in de Sitter space. 
This state is de Sitter invariant, but each geodesic observer sees a thermal bath in her 
rest frame with the same universal temperature T = H/ (2tt) where H is the inverse of the 
curvature radius of the de Sitter space. 

Quantitatively, we find that the minimum uncertainty in the determination of the time 
of nucleation At nuc , which can be inferred from the response of the detectors, is of the 
order of the size of the instanton (see Eq. (7.16)). It therefore appears that this can be 
much smaller than the timescale Ai m i X i ng for mode mixing (estimated from Eq. (2.25)). 
For A = m^/(eE) > 1 we have At nuc ~ A^ 1/2 At mixing . 

s This scenario was suggested by N. Tanji [14], although he did not provide evidence to support it. 
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The close analogy between pair creation in (1 + 1)D and bubble nucleation leads us to 
expect that a similar picture should apply to observation of bubble nucleation (in (3+1) as 
well as in other dimensions). An observer O who fills the space with detectors at rest with 
respect to himself will see bubbles nucleating at minimal radius at rest and then expanding 
in the outward direction. Another observer O' , who has installed detectors at rest in her 
frame, sees the same picture. But if O' were to watch the detectors of O, she would see 
a very different scenario. A small piece of the bubble wall first appears and moves at a 
very high speed towards the bubble center. It gradually slows down and grows in size, and 
eventually bounces at the minimum radius when the wall extends over half a sphere. After 
the bounce, the wall expands, and still later it closes up to a full sphere. 

We also discussed the unusual, and to a certain extent unphysical, properties of the in- 
vacuum state, which put the above stated conclusions somewhat in doubt. The singularity 
structure of the two-point function in the in-vacuum state does not have Hadamard form, 
and as a result the expectation values of physical observables cannot be regulated in a 
Lorentz invariant way. The same could apply to regularization of the loop diagrams in 
higher orders of perturbation theory. This violation of Lorentz invariance has no effect on 
the tree-level scattering processes that we discussed in this paper. However, it does raise 
some questions about the overall consistency of the model. 

In a more realistic model, the electric field would have to be switched on at some initial 
time in the past. In order to investigate the effect of initial conditions, one would have to 
do similar calculations for an electric field which is turned on for a finite period of time T 
and check whether or not the results approach our results in this paper in the limit T — > oo. 
It would also be interesting to extend the analysis to de Sitter space 9 . We hope to address 
these problems in future work. 
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A Klein-Gordon product 

Here we calculate the Klein-Gordon product of mode functions, defined in Eq. (3.15). Using 
the integral representation (2.15) for the parabolic cylinder functions, we can write, up to 

9 Schwinger pair production in (l+l)-dimensional de Sitter space has been considered in Ref. [22]. 
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normalization, 



2 fOO 2 

<j) k {t,x) = e iS * toe^-T /, (A.l) 
J o 

and 

-2 Z* 00 /2 

fa(t,x) = e- ien e iS 2 dw' ea-iW-V (A.2) 
Substituting in (3.15) we have 

/•OC 2 /"°0 ,2 

1C= / dww u e~^ / dw' w /u e~^I(k,k,w,w',t), (A.3) 
Jo Jo 

where 

J(fc, fc, w, w',t) = i (eE) 1 / 2 [ dx [(1 - vhz -z-(l + i)( 7 w' - w)]fj- k (A.4) 



J t=const 

and 

f k = e ~ ikx e i2 ^ e {1 ~ i)zw , f~ k = ^-ieO e (i-*)«e' _ ( A _ 5 ) 

Let us consider the mixing of the positive with the negative frequency modes. For 
convenience, we choose the integration surface at t = and take k = 0. In this case, 
Eqs. (3.14) give x = jx, z = —^feE^vx + k/yeE, and we have 

I = _L e *35SB f dz[(l-vhz-(l+i)hw' -ujje-^^N™' . (A.6) 
7^ A=o 

Performing the Gaussian integral, we have 

2\/iri i k 2 vw' 2 

1= — v e 2 ^ E (w - 7(1 + v)w') . (A.7) 

7V (1 - v ) v 

Finally, performing the w and w' integrals, we obtain 

K = . (A.8) 

B Two-point function 

In this Appendix we shall calculate the gauge-invariant 2-point functions G^, defined by 
Eqs. (4.5), (4.6). We consider the in- vacuum expectation value, 

G + (zV) = in <O|0V) e-feJ?^" 0(^)|O) in . (B.l) 
Using the expression for the 0-field given by Eq. (2.8), we find 



iefA„dx» / dk ^ e ik(y-x) _ ^ ^ 
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Using the integral representation of cftk from (2.15), the integral over momenta leads to a 
delta function, and we have 



G+ = ]^! e --J-^ I dX X^e^ 2 I " dX' X'^e- 1 ^' 2 
4vr J J 

xei eE(ti-t>i)-iV^E(tx-t'x>) 6 (x'-X- 5_) , (B.3) 
where a is the Bogoliubov coefficient given by Eq. (2.22). We also define 

FeE 

S± = J — (At ± Ax) , At = t' - 1 , Ax = y-x (B.4) 
Assuming 5_ > 0, we have X' > X. Doing the X' integral we obtain 

I 1 2 r 

4tt J 

(B.5) 

This can be rewritten as 

G + = rfj-fe/M'^lC-O+i^^-^. J dx{x + S jfx» e iS + x > (B _ 6) 

and, after some rearrangements, 

G+ = ^e-^S A v dx» e -^eE&x{t+t>)^_ iS+s _y\ T{-u*) W iX/2 (-iS+S-) , (B.7) 

where W a> p (z) is the Whittaker function. 

In the gauge (2.6), if we integrate along a straight line, so that dx = (Ax/ At)dt, then 
the Wilson line is 

/ A v dx v = -—Ax(t + t'), (B.8) 
so that the two-point function is given by: 

lal 2 W iX/2 n (ieEAs 2 /2) 

G + (Ax") = l -^-F(-u*) lX/2 ' ° 1 1 L , for At - Ax > (B.9) 

4vr y/ieEAs 2 /2 

where As 2 = -(At) 2 + (Ax) 2 . 

A similar calculation for S_ < gives the complex conjugate expression: 

Id 2 W-ix/2 o He£As 2 /2) 

G+(Ax^) = ^r(-i/) tA/ 2 '° l — L^-, for At-Ax<0 (B.10) 

47r y/-ieE As 2 /2 

Likewise, we may consider the 2-point function 

G" K, y") = in (0| cf>(yn e~ ie f* A » dx " V) | ) in , (B.ll) 
which can be expressed in terms of the mode functions as 

2vrG"(At, Ax) = e -™J A ^» J dk ^{t') (^(t) e ik ^~ x) . (B.12) 

Comparing this with Eq. (B.l), we see that it is simply related to G + , 

G-(At,Ax) = G + {-At,Ax). (B.13) 
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C Lorentz violating current 



In this Appendix, we evaluate the expectation value of the current in the in-vacuum, 
using a direct momentum space regularization. 

Substituting the mode expansion (2.8) in Eq. (4.3), we have (Jo) = and 

(Ji) = 2e J ^(k + eEt)\4> k \ 2 . (C.l) 

Here, <j) k are functions of z = {eE)~ l / 2 [k + eEt), so the integrand in Eq. (C.l) is a function 
only of z. It appears that we can perform a change of variable k —> k — eEt to make 
the integral independent of t. However, such a shift of integration variable is not always 
legitimate for divergent integrals. 
Let us consider 

This shows that (dji/dt) is nonzero if the function 

F(k,t) = (k + eEt)\(j) k \ 2 (C.3) 

has different limits at k — > ±oo. 

Now, the asymptotic forms of 4>k(z) are given by 

4>k(z) « f k {z) , for z ->• -oo (C.4) 

and 

<f> k (z)naf* k (z)+Pf k {z), for z -»• -oo, (C.5) 

where 



f k (z) « (e£)-3 -^r-r- (C.6) 



(v^kl) 

Substituting this in (C.3), we have 

F (k -> -oo) = i (C.7) 

For fc — >■ +oo, -F(/c) does not approach a definite limit, as it includes rapidly oscillating 
terms: 

F(k -)■ +oo) = - (|a| 2 + |/3| 2 ) + (oscillating terms). (C.8) 

The oscillating terms indicate that dJ\/dt is not a well defined quantity in the in- 
vacuum. In order to remove these terms, let us consider 

(Ji(t) - Ji(O) = jf* «*r^(r) = ^ jf [(k + eEr)h fe | 2 ]^ . (C.9) 

We can evaluate this by first integrating over r and then taking the limits k — > ±oo. This 
corresponds to the usual subtraction procedure, where we first do the subtraction in the 
integrand and then integrate over momenta. 
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The oscillating term is proportional to 

(CIO) 



» 2 2ifct 



This oscillates infinitely fast as k — > oo, so the oscillating term in (E.7) vanishes in the 
limit, and we obtain 

(Ji(t) - Mt')) = —e~^(t - t>). (C.ll) 

Here, we have used 

|a| 2 = l + |/3| 2 (C.12) 

and 

|/3| 2 = e-^. (C.13) 
Finally, taking the limit t' — > t, we obtain 

dJA e 2 E _™£ 

■w) = — 6 eE - (ai4) 

D Transition amplitude for a ^-detector 

Here we shall calculate the integral over mode functions appearing in Eq. (5.26) for A x (p; q) 
following the method developed in Ref. [13]. 

We first express the mode functions (2.14) for <p and ip fields using the integral repre- 
sentation (2.15) of parabolic cylinder functions. The neutral x-P ar ticl e mode function is 
simply given by a plane wave. Thus we have 



e 2 <t> 



dX X"* e~ i^ z 4> X + 







™ (27r) 1 /2(2e J B)i/4 

* - wvw* ^syy* 

where the Bogoliubov coefficient a defined in (2.22) is used with an index of each particle 
instead of its momentum, because we found that there was no dependence on momentum. 
We have also defined 



eE[t + - , = _, K ^ 



^E[t + —) , ff = , A ^ = ^' ( D ' 2 ) 

Plugging Eqs.(D.l) and (D.2) into (5.26), the amplitude A x (p;q) becomes 



A (va) - ga<t> ^ e^iP-P 2 ) / dX X u l Jv ^ {q ~ p)X ~ i: 



' J pr jo 

,2 -I POD 



pOO /2 V 2 1 POO 

x / dY e - i y^ qY+l - — / dt e i(p+^+v^E{x-Y))t _ (D 3) 

J0 2vr J-oo 
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The time integration gives rise to the delta function 

5(V2eE(X -Y)+p + 



(D.4) 



The delta function tells us that Y > X, and so we can perform the Y integration first. 
Then, we have 



where 



2eE (2w p ) 1 /2 



~ I dX X v * (X + fi + ) v * e ifl - x (p.5) 
o 



Up ±p 



(D.6) 



The remaining integral can be expressed in terms of the Whittaker function. By using 



m x _ A x 



2eE 2 



we finally obtain 

l-4x(p;?)l 

where 



tt \ 1/2 g\foa%\ 



Up J \/2eE 



m v 



Wur +t ia_ | 1 , 



•V ± A,/, 



(D.7) 



(D. 



(D.9) 



and 1/3^1 = e ^Z 2 is the negative frequency Bogoliubov coefficient. 

E Transition amplitude for a %-detector 

Now we calculate the transition amplitude A^(q;p) in (5.28). 
Using the identity 



D v * [-{l-i)z\ 



T(is* + 1) 



e^ v LU*_i[-(l + i)z] + e"3 w + i)z] (E.l) 

and the integral representations (D.l), we have, up to an overall phase factor, 

/oo 
dt 4>* q _ p ipl x P 
-oo 



9a % 



(2e J B) 1 /2(2 Wp )i/2 



dY Y v *if e J V eE 



i7 ™<t> e 2h(p - 2 ip) I dX X^ e~V ^ 



{q-p)X + i 



2 



2tt 



dt e i(-p-u> p -V2eE(X-Y))t 



+ 



90% 



(2e£?)i/2(2o; p )V2 7 Q 



: / dY Y v + e V 
/o 



2?r 



/oo 
dt e i{-p-uj p +VteE(X+Y))t _ ^ E _ 2 ) 
-oo 
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The time integration in the first term in Eq. (E.2) give a delta function 

I — p — UJp — V2eE{X-Y)^ . (E.3) 

This delta function tells us that Y > X, so we can perform the Y integration first. Then 
the first term is 



2eE (2w p ) 1 /2 
where 



. poo 

2 



rfir* (x + e -^- x , (e.4) 

(E.5) 



On the other hand, the time integration in the second term of (E.2) gives 



5 l-p -uj p + V^E(X + Y)j . 
Performing the Y- integration first, we find 



3«V 



2eE(2u p ) 1 / 2 



dX X"+ (-X + n+y* e 



< Jn-x 



(E.6) 



(E.7) 



The remaining integrals of Eqs.(E.4) and (E.7) can be expressed in terms of Whittaker 
functions. Using 



fjt+fi. 



m x — 



2eE 



(Ei 



we finally obtain 

\Aj>(q;p)\ 



TT 



OJ, 



1/2 



a. 



y/2^E 



W, 



ICT+, 1(7- 



4 1 



where 



A</> ± A.0 



2CT+, ?tr_ 



(E.9) 



(E.10) 



F Amplitudes for a time-dependent coupling 

To evaluate the amplitude A x (p;q = 0) in Eq. (7.8), we first substitute the plane wave 
solutions for ipo and Xp an d the integral representation for 



/oo 
dt g(t) <P-p(t)Mt)x* P (t) 
-oo 



Cu~ 1 / 2 I dXX v *e 



2 /*OC 



dte t 7 e * 2 e * 



V2zX 



where 



(F.l) 
(F.2) 
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and 

v ^F^(2e^) 1 /4 1 ' ^ 

is a constant. After completing the square in the exponent of the i-integral and performing 
the Gaussian integration, this becomes 

A x (p; q = 0) = CT^ 1 ' 2 ^ - ^ ET2 e ™ "^T jf dXX v *e^ x ^ X , (FA) 
where 

= 2^+ y/2e£r 2 (^-m^) 

The integral in the last equation can be expressed in terms of a parabolic cylinder function 
using Eq. (3.462.1) of Ref. [21], 

jf dXjtf e-W-i* = V{1 + /^-^D.^ ^ J . (F.7) 

We are interested in the p-dependence of the amplitude. From Eqs. (7.2) and (7.7), it 
is easy to show that the regime of interest to us is \ r y/y/2j3\ 3> 1, and using the asymptotic 
form of the parabolic cylinder function we have 

e z ' 2/i D v {Z) oc Z v . (F.8) 

Then it is easy to see that the dominant factor determining the p-dependence in this regime 
is the exponential factor in Eq. (F.4), 



exp 

Hence, 



T 2 (p + uj p -m^ 



\A x (p;q = 0)| 2 cx exp 



2(2 + ieET 2 ) 

(p + uj p - m^) 21 



(eETf 



(F.9) 



(F.10) 
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